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We show that how a recent experiment of quantum imaging with undetected photons can basically be de¬ 
scribed as a (partial) ancilla-assisted process tomography. We propose a simplified quantum circuit version of 
this scenario, which also enables to recast quantum imaging in quantum computation language. Our analogy and 
analysis may help better understand the role of classical and/or quantum correlations in imaging experiments. 
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I. INTRODUCTION 

Quantum imaging (QI) is an interesting technique which 
employs quantum properties of light, such as entanglement 
and nonlocality, to give a high resolution image from a 
partially-transmitting object [1]. QI has also been experimen¬ 
tally demonstrated in numerous experiments [2], e.g., through 
the method of “ghost imaging” [3-6]. 

To recover information about an unknown ‘object’ (illu¬ 
minated by either a “signal” or an “idler” photon), most QI 
methods employ entanglement of quantum states generated 
by parametric down-conversion. However, it has been shown 
that entanglement is not necessary for ghost imaging since us¬ 
ing classically-correlated fields [7], thermal [8], and pseudo- 
thermal lights [9, 10] still enables obtaining images through 
ghost imaging. In addition, it has been shown that using half¬ 
wave plates properly (considering polarization as a degree of 
freedom), one can affect output probabilities and thus visibil¬ 
ity of the images [11]. Recently, a novel and elegant method 
for QI has been experimentally demonstrated in Ref. [12], 
where “which-path entanglement” seems to be crucial in or¬ 
der to have image. A closer theoretical inspection of this ex¬ 
periment from the perspective of quantum optics has been re¬ 
ported in Ref. [13]. 

Quantum process tomography (QPT) is another technique 
employed to identify unknown quantum “processes” [14]. 
Several schemes have been outlined to accomplish this task, 
namely, standard QPT [14-16], ancilla-assisted process to¬ 
mography (AAPT) [17-20], and direct characterization of 
quantum dynamics [21-23] —see Ref. [24] for an extensive 
review. Notwithstanding their differences, all QPT methods 
operate similarly based on probing an unknown process (as a 
“black box”) with appropriate input states, measuring output 
states, and identifying the process through relation of input 
and output states. 

Establishing a connection between quantum tomography 
and imaging or sensing schemes can be interesting from var¬ 
ious aspects. It has been argued that tomography and spec¬ 
troscopy can be considered as dual forms of quantum compu¬ 
tation [25]. Following a similar reasoning, one may also argue 
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that QPT and QI of an object are basically akin in the sense 
that they both examine a(n) process/object with probe states 
and then analyze output states. Here we make this connec¬ 
tion more explicit. In particular, we demonstrate that the QI 
scheme proposed in Ref. [12] can be described as a version of 
AAPT in which the object is assumed in both as a black box, 
whereas an identical or partial “image” is obtained by ana¬ 
lyzing ancillary probes. We elaborate in detail how this anal¬ 
ogy work through identifying different steps of the QI scheme 
with preparation and measurement parts of a special AAPT 
scenario. Additionally, we represent a quantum circuit for the 
QI scheme, which helps analyze the role of initial-state corre¬ 
lations in QI. We demonstrate that, by replacing the entangled 
Bell states with arbitrary Werner states, QI can also work by 
non-entangled light. 

This manuscript is structured as follows. In Sec. II, we re¬ 
view the QI scheme of Ref. [12]. In Sec. Ill, we suggest a 
quantum circuit version of the QI scenario in which only one- 
qubit and CNOT gates as well as a measurement are used. In 
Sec. IV, we describe an analogy between the QI scheme and 
AAPT by translating steps of the QI scheme into a version of 
AAPT in which parts of the probe systems remain undetected. 
We explicitly show that which parts of the QI scheme cor¬ 
respond to preparation and measurement parts of the AAPT 
scheme. In addition, we show that quantum entanglement of 
probe states is not essential in generating an image, and that 
how one can supplement measurements to enable a full to¬ 
mography of the object. We summarize our findings in Sec. V. 
There are several appendices which include details of parts of 
calculations. 


II. THE IMAGING SETUP 

Here we focus on a recent method and experiment for cre¬ 
ating “image” of an object in a quantum mechanical fashion, 
proposed in Ref. [12] and theoretically furthered in Ref. [13]). 
Our discussion here concerns (a simplified version of) this 
experiment without realization technicalities, and aims at ex¬ 
plaining how imaging an “object” can be understood in the 
framework of QPT. 

In this experimental setup a beam splitter (BS1), which is 
illuminated by a pump photon, is used to generate the follow- 
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FIG. 1. (color online). Experimental setup of the QI scenario in 
Ref. [12]. 


ing path-entangled Bell state: 

l^l) = (|l>p,a|0) p ,b + |0> p ,a|l) Plb )/^, (1) 

where p indicates the wavelength of the pump photon, and 
a and b denote different paths which the photon can choose 
(Fig. 1). To set our notation hereafter, we assume that | n) xt 
indicates an n-photon state with wavelength A in path t. Since 
at any instant there exists at most one photon in each mode , we 
use the encoding \n = 0) = (1 0) T and \n = 1) = (0 1) T , 
which represent a logical basis as a “qubit” [14]. 

This setup includes two nonlinear crystals used to generate 
extra photons. These crystals should be pumped identically 
but not simultaneously. While passing through the nonlinear 
medium (NL) [26], a pump photon (532 nm, in the original 
experiment) can convert to a signal (s) and an idler (i) pho¬ 
ton with different wavelengths (810 nm and 1550 nm, respec¬ 
tively, in the original experiment) through a down-conversion 
process. This is a probabilistic event. Hence the probability 
of obtaining a photon in the final detectors of the setup are 
conditioned on a successful down-conversion process. How¬ 
ever, in the following, we do not consider unsuccessful down- 
conversion events as they do not yield an image in this setup. 
NL1 changes the Bell state (1) to 

|*2> = i(|l) 8 ,a|l)i.a|0)p,a|0) p>b + |0> i , a |0) i , a |0> p , a |l) p , b ). 

( 2 ) 

The signal and idler photons are next separated (because of 
their distinct wavelengths) by a dichroic mirror (Dl), which 
reflects the idler photons into path d and allows photons with 
other wavelengths to pass. This yields 

1^3) = y^(|l)s,c|l)i 5 d|0)p,c|0) Pj b + |0) SjC |0)i J d|0)p >c |l)p 5 b)- 

(3) 


of Ref. [12], the object is placed between two lenses, and it has 
been argued that [13] an incident plane wave mode with wave 
vector k can illuminate one point of the object and hence 

yields an output wave with an associated wave vector k' ( k). 
Given that a quantized light is a superposition of plane wave 
modes, one can assume that one point on the object can trans¬ 
mit and reflect only one specific mode of the quantized (idler) 
beam [13]. Bearing in mind this one-to-one correspondence 
between the object and the image, it thus suffices—for our 
purposes here—to explain the basic theory of the imaging sce¬ 
nario through a single point of the object with constant T and 
7 factors. One may consider this restricted sort of imaging as 
a “sensing” scenario. 

After interacting with the object, |\h 3 ) changes to 

|*4> = A; (Te i7 |l) S)C |l) i)d |0) i)W |0) P)b |0) p , c 

+ Vi - T^|l).,c|0)i,d|l)i,w|0) p ,b|0) p , c 

+ |0) s?c |0)^d|0}i,w|l)p,b|0}p, c ), (4) 

where w is the path through which the reflected idler photon 
passes. There is a probability that after BS1 the photon goes 
through path b and the whole state be affected by the dichroic 
mirror D2. This mirror maps paths d, b —> e, thus the state of 
the system is transformed to 

A ; (r e < T f |l) i , c |l) i ,e|0) ilW |0> Pl e|0) p , c 

+ x/l - r 2 '|l) S) c|0) i) e|l) i ,w|0) p ,e|0) P)C 

+ |0) s , c |0)i ) e|0)j iW |l) J , i e|0)p, c ) . (5) 

Now NL2 acts on this state as |0) Sj e|0)i, e |l)p,e —>* 
|l) S5 e|l)i,e|0)p, e ; that is, it transforms a photon with the pump 
frequency to two photons along the same path e. Hence the 
total state becomes 

1*5) (Te < T r |l) i , c |0) 8) e|l>i,e|0)i )W |0) p ,e|0) p , c 

+ 7l- T 2 '|l) 5 , c |0) 5 ,e|0) i) e|l) i> w|0) P) e|0) P)C 

+ |0)s,c|l}s,e|l)j,e|0)j > w|0)p ) e|0)p jC ) • (6) 

By using another dichroic mirror (D3), the idler photons are 
discarded from the setup. The pump photons are also dis¬ 
carded through D4 and D5; the reflected photon from the ob¬ 
ject is discarded too. Mathematically, such discarding is rep¬ 
resented by tracing out the idler photon, the pump photon, and 
the photon reflected to path w. Hence, the total state before 
BS2 becomes 


The idler photon then illuminates the object located in path 
d; it either passes through the object (O) with the associated 
transmission factor Te* 7 , or is reflected with the associated 
reflectivity factor y/l — T 2 . 

We note that for a real two-dimensional object both T and 
7 depend on position on the object, that is, T = T(x, y) and 
7 = 7 (x,y), where (x,y) G O. In the real experimental setup 


26 =|(|l)c<l| ® 10) e (01 + Te*\ l)e(0| 0 10) e (11 

+ Te-70) c (l| 0 11) e (01 + 10) c (01 0 |l)e(l|), (7) 

where to ease the notation we have removed the wavelength 
indices because after tracing out the idler and pump photons 
only signal photons remain. 
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A pivotal feature of this setup is that after NL2 the source of 
the idler photons cannot be distinguished, and this fundamen¬ 
tal ambiguity in “which-path information” is indeed respon¬ 
sible for creating the image. In other words, |l)z,d = |l)i, e » 
which indicates that after NL2 the source as well as the path 
of the idler photons cannot be discerned. At the final stage, 
BS2 is applied to combine the information carried in paths c 
and e, and to generate the image. The probability of finding a 
photon in paths h and g is obtained as 

Ph/ g = (lTTcos 7 )/2. (8) 


III. THE IMAGING SETUP REPRESENTED BY A 
QUANTUM CIRCUIT 


In this section, we present a quantum circuit which can sim¬ 
ulate the setup of Fig. 1. 

We assume a quantum circuit as in Fig. 2, which is ap¬ 
plied on four qubits prepared in the different field modes 
|$o> = | 0 )si| 0 >ii| 0 )i 2 | 0 > S 2 , where indices “1” and “2” imply 
the source of the qubits (or photons), namely, NL1 and NL2. 
The first Hadamard and the three subsequent CNOT gates turn 
|$o) to (App. A) 


l^probe) = (|l) 8 l|l)il| 0 ) i 2 | 0 ) a2 + | 0 ) 8 l| 0 ) j i|l)i 2 |l> a 2 )/v^. 

(9) 

The Hadamard followed by the first CNOT (which acts when 
the control qubit is |0)) play the role of BS1, and the other two 
CNOTs simulate the action of NL1 and NL2 in the experimen¬ 
tal setup. 

The object (O) is represented by the quantum channel (pro¬ 
cess or operation) Eq which acts as (App. B) 


rio)(oi^io)(oi 

I |0}<l|^Te- 7 0)<l| 

°‘ | |1)(0| -> Te i7 |l)(0| 

[|1)(1|^T 2 |1)(1| + (1-T 2 )|0><0|. 


( 10 ) 


After applying we need some operation which can make 
|l)i,d = |l)z, e in the original experimental setup. In our cir¬ 
cuit, this can be achieved, for example, if we make the two 
states 11)?! 10) 2 2 and |0)i 1 11)» 2 quantum-mechanically indis¬ 
tinguishable. To do so, we can consider a mode mixer (MM) 
such that 

MM: <n) 




|0)a - 

l°)*2 - \ji\- 

|0> a2 - 




—i 

u 




D 




-o 


€^[> 


FIG. 2. Quantum logical circuit of the imaging setup of Fig. 1. Here 
H is the Hadamard gate and Z is the z-Pauli gate. 


and acts as identity otherwise, where |H) is a fixed state of the 
idlers (App. C). The MM can be realized in various manners. 
As an example, consider one Hadamard gate, one cz gate, 
and the projection n + = |+)(+|, applied respectively, where 
| =b) = (1/^/2")(10) zb 11)). For this specific operation we have 
|S)iii 2 = |—)iil+)i 2 - I n general, the action of MM can be 
described through the following quantum operation: 

«r r,_ (PH + Q)P(PH + Q) f 

SmmIq\ Tr[(P=; + Q)^(Ph + Q)t] ’ ( ’ 

where P- = |S)((01| + (10|) and Q = (|00) + |11))((00| + 

(ni)- 

The action of MM and tracing out over the idler photons 
(simulating their absorption or loss in Fig. 1) reduce |<f>p ro be) 
to 

T =i(|l)si(l| 0 |0) s2 (0| + Te* 7 |l) sl (0| 0 |0} s2 (l| 

+ Te _ * 7 |0) sl (l| 0 |l) s2 (0| + |0) sl (0| 0 |l) s2 (l|), 

(13) 

which is equivalent to the state qq [Eq. (7)] in the original 
experimental setup. The details can be found in App. C. 

The action of BS2 can be simulated by a Hadamard and two 
CNOT gates on T. To detect photons in the first and fourth 
paths with correct probabilities P Sl /s 2 = (1 T Tcos 7)/2 
[Eq. (8)], we shall need an appropriate measurement fW, 
which we discuss in the next section. 


IV. QI VS. AAPT 

Here we revisit the QI scenario and demonstrate that it 
can be captured as a form of AAPT. We first explain a gen¬ 
eral framework for AAPT with undetected particles. Next we 
translate the QI setup into this modified AAPT language. Here 
by “undetected” we mean that we discard part of the whole 
system and employ measurement results on the remained parts 
to infer the unknown channel. 


A. AAPT with undetected particles 

In the standard picture of an AAPT scheme [17-20], an 
unknown quantum process (or black box) acts on a specific 
system which is correlated with an ancillary system. Iden¬ 
tification of the process comes through analyzing results of 
measurements on the whole composite systems [24] . 

Another variant of AAPT may include undetected sys¬ 
tem^), where measurements are instead performed only on 
ancilla(s)—whereas the rest of the systems are discarded. In 
such indirect schemes, although measurements are performed 
on systems that never directly met (i.e., passed through) the 
black box, because of correlations between the system and 
ancilla(s), information of the black box would still leak to the 
ancillas such that one can identify the black box by measur¬ 
ing the ancilla(s). In the following we briefly describe this 
modified version of AAPT. 
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FIG. 3. An AAPT scheme with unmeasured particles (discarded). 
Here 9V[ denotes measurement on two of the ancillas. 

An arbitrary state of a bipartite system (, system + ancilla), 
with the Hilbert space !Hi 0 (as input of AAPT) can be 
written as (Schmidt decomposition) 

D 2 

g (m) = rt A t (g, (14) 

1 =1 

where m s are some (nonnegative) numbers, and the 
({^}|Li) f° rm an orthonormal operator basis on 9^ (9( 2 ), 
satisfying the orthonormalization condition Tr^Aj,] = <5^/ 
(TrfF^Dj,] = Suf), with D being the dimension of the small¬ 
est Hilbert space [27]. It suffices for tomography to choose 
the dimension of the ancilla the same as the dimension of the 
system. After the black box £ acts on the system, the output 
state of the total system becomes 

2 (out) = (£ ®X)[2 (in) ] = ® B e , (15) 

£ 

where 1 denotes the identity operation. In the standard AAPT, 
measurements on p( out ) yield a set of linear equations by solv¬ 
ing which one can read £. We now modify this picture as 
follows. 

The state of the ancilla can be read from p( out ) by tracing 
out over the state of the system, 

02° Ut) =XV/'fr[!£[i4 / ]]B / . (16) 

£ 

This yields 

{Bl) out =r e Ti[‘Z[A e }], (17) 

where (s]) ou t = Tr[B ( t £> 2 ° ut ^. If Be s are Hermitian oper- 
ators, this relation gives the result of measuring the observ¬ 
ables Bi on the final state of the ancilla. If the number of 
nonzero 77 s (i.e., the Schmidt rank) of the initial state is equal 
to D 2 and the number of independent parameters of £ is not 
greater than D 2 , this system of linear equations can be solved 
[19, 24]. Thus in this particular case, we see that measurement 
on the ancilla would be adequate to reproduce the unknown 
process £. 

We note, however, that a general quantum chan¬ 
nel/operation acting on a system with a D-dimensional Hilbert 
space has D A — D 2 independent parameters [14]. Hence the 
above modified picture does not apply to general processes; 


A B 



FIG. 4. (color online). A circuit-like version of the QI process of 
Fig. 1 . Here parts A and B show the state preparation and measure¬ 
ment parts of the corresponding AAPT. Note that Fig. 1 does not in¬ 
clude the phase shifter 0 (in the measurement part). The role of this 
extra device is to enable full tomography of the object operation— 
see Subsec. IV C. 


it applies only when the number of independent parameters 
in £ is not greater than D 2 . An example of such restricted 
case is depicted in Fig. 3. Here we have four systems; one 
on which the (unknown) process £^ acts, and three extra (an¬ 
cillary) systems, all with the same Hilbert space dimension 
d. For this setting the process £ of Eq. (15) is in the form 
£ = £0 which has d 4 — d 2 parameters. That is, here we 
have D = d 2 . Note that in this modified AAPT we discard 
the system (one on which £o acts) and one of the ancillas but 
measure on the remaining two ancillas. In Fig. 3, 9/[ denotes 
measuring the {B#} observables. 

As a remark, note that if the Schmidt rank of the initial 
state is R (obviously ^ D 2 ) and at the same time the number 
of unknown parameters of £ (exactly speaking, £q) is ^ R , 
the modified AAPT method still works. 

The above modified AAPT can be generalized further. One 
can consider that before discarding the undetected parts, a 
known operation is applied on them. Since a quantum op¬ 
eration is linear, the above argument about solving Eq. (17) 
still holds but now the equation is modified as 

(Bl) out /r i = Tr[fo‘E[A e }], (18) 

where “o” denotes composition of quantum operations. 


B. QI as an AAPT 

A circuit-like version of the QI setup is sketched in Fig. 4, 
where the dashed parts A and B represent, respectively, the 
state preparation and the measurement [24] . 

For the imaging setup of Fig. 2 (or Fig. 4), £o (O) acts 
on ii (id) and has only two independent parameters T and 7 
(App. B), and the idlers are undetected (hence D = 2 2 ); we 
only detect the signal particles. Here the initial state is the 
entangled state (9), 

4 

@iii2,s 1 S2 = l^probe) (^probel = ® (^) SlS2 , 

£=1 

( 19 ) 
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= r 3 = — r 4 = 1/2 and 


= Pi = (1/V8)(II-ZZ), 

(20a) 

= B 2 = {l/V8)(zi - IZ), 

(20b) 

= Bs = (l/V8)(XX + YY), 

(20c) 

= B 4 = (1/Vs)(XY - YX). 

(20d) 


Here for brevity we have removed the tensor product symbol 
(thus, e.g., XY = X 0 Y), and X , F, and F are the Pauli 
matrices [14], and I is the identity operator. Equation (19) 
is the Schmidt decomposition of g 6 n ). But note that here we 
have chosen r 4 = —1/2 in order to make { A 1 } and {B 
(orthonormal) Hermitian operators. Although p( m ) is not full- 
rank (R = 4 < 2 4 ), it suffices to determine £o- 

After applying the object operation £0 and the MM oper¬ 
ation 5mm on i ii 2 , the state of the total iii 2 siS 2 [Eq. (16)] 
becomes 


S =2H 


?mm ° (£0 0 !)[// - ZZ] 0 (II - ZZ) + ?mm o (£0 ®I)[Z/ - /Z] 0 (ZJ - JZ) 

+ Jmm o (‘Eo®i)[xx + yy] 0 (xx + yy) - j M m ° (£0 0i)[xy - yx\ 0 (iy - yx) 

V2 


=Y([(1 + r 2 )l=>(2| + (1 - t 2 )| 00 ) 00 |] 0 ^ + [(1 - T 2 )|S)(S| - (1 - r 2 )|oo)oo|] 

, w . XY-YX , w , XX + YY\ 

- 2 T sin7 |=.)(=.| 0 --+ 2 T cos 7 |s)(s| 0 --J. 

Discarding the i\i 2 systems (i.e., tracing out over iii 2 ) yields the state of the signal systems (5182) as 

XX + YY m m XY-YX 


ZI-IZ 


T = 


X 

2 


II-ZZ 

~W 


+ T cos 7 ; 




Tsin7 ; 




( 21 ) 


( 22 ) 


This relation indicates that in order to obtain the object pa¬ 
rameters T and 7, it suffices to measure the observables 
XX + YY and XY — YX. These measurements conclude 
full AAPT of the object. In the quantum circuit of Fig. 2, the 
last part (a CNOT, a Hadamard, and another cnot) before the 
detectors, represent BS2 in the experimental setup of Fig. 1 or 
the beam splitter in part B of Fig 4. The operators 


M h = |l) h <l|0|O) g (O|, 

= |0)h(0| (8 |l) g (115 


(23) 

(24) 


indicate the action of the final two detectors—that is, we ei¬ 
ther detect a photon in path s 1 (h) or a photon in path 82 (g). 
Hence, the measurement performed on 8182 (hg) can be de¬ 
scribed by the following two observables: 


fTf h/g = (cnot H <g) / CNOT^ M h /g (cnot H ® I CNOT^ 

(25) 


II-ZZ XX + YY ] 
T 


5s Vs 

It is straightforward to see that indeed are the following 
Bell-state measurements (App. D): 

A4/ g = I^X^I, (26) 

where l^) = (|01) =F |10))/v^. 

With these Bell-state measurements, the probability of de¬ 
tecting a photon in 81 (path h) is obtained as 


P h = TV[fV h T] = (1 - Tcos 7 )/2, 


(27) 


and similarly for detecting a photon in 82 (path g), 

P g = Tr[5V g T] = (1 + Tcos 7 )/2, (28) 

in agreement with Eq. (8). 

Several remarks are in order here, (i) It is seen that 
Ph + P g = 1, that is we definitely detect one photon either 
in path h or in path g. (ii) It is evident that with only Ph or 
P g (either or even both), or equivalently the Bell-state mea¬ 
surements one cannot obtain the complete information 

of the object (i.e., both T and 7 at the same time). Hence 
the QI setup of Fig. 1 corresponds to a partial AAPT. As we 
commented after Eq. (22), one further needs the other (non- 
Bell-state) measurement XY — YX to fully characterize the 
object O. We discuss later how one can realize this measure¬ 
ment by a simple modification of the original setup of Fig. 1. 
(iii) It is clear that the initial probe state |4> pro b e ) [Eq. (9)] is 
an entangled state. We also discuss below how essential this 
entanglement is for the QI setup to work. 


C. Complete tomography of the object 

Here we focus on the dashed part B in Fig. 4, which repre¬ 
sents measurement in the QI setup. As we discussed earlier in 
this section, the final measurement in the QI setup of Fig. 1, 
which gives the probabilities P h / g , is tantamount to a partial 
tomography of the object—in the sense that through this mea¬ 
surement we only obtain the value of T cos 7, not both T and 
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FIG. 5. Probability of detecting a photon in path h as a function of 
the phase parameter £>. This plots shows that how by varying £> one 
can read both T and 7. 


7 separately. We also argued that a full tomography would 
require another measurement (XY — YX). Here we demon¬ 
strate that in order for the full tomography of the object, one 
needs to slightly modify the existing measurement by adding 
a controllable phase shifter (£> in part B of Fig. 4). We show 
that this modification enables measurement of the missing ob¬ 
servable. 

One can see that the action of the phase shifter is to change 
fWh/g to the following observables (App. E): 


r ((k\ V2 r n — zz ( 

, XY-YX\ 


XX+ YY 


(29) 


Hence, in the output of the detectors we obtain the following 
probabilities (rather than P h / g ): 


Pi% = Tr T] = (1/2) [ITT cos( 7 + 4>)\ ■ (30) 


It is straightforward to see that by varying £>, one can read 
both T and 7 through (or Pg^) as in Fig. 5, which is a 
technique used in the Mach-Zehnder interferometry [28]. 


D. Role of quantum correlation in QI 

It may seem that quantum entanglement in the QI setup of 
Fig. 1 is essential to obtain an image. However, employing 
the AAPT version of the QI setup allows to see that even sep¬ 
arable initial states can be useful in order to extract an image. 
The fact that entanglement is not essential has already been 
proved for AAPT [19]. 

As seen from Fig. 2, the initial (prepared) state right before 
the object is |4> pro be) [Eq. (9)], which is an entangled Bell 
state, which can be represented in a form more familiar for 
the Bell states as 


| ^probe) 




+ |l)si*l |0)»2S2)> 


(31) 


by using the encoding |0) Smim = |0) s j0) im and |l) Smim = 
(me {1,2}). 

We replace this (extended) Bell state with the following (ex¬ 
tended) Werner state [29] 

W"ext = (C/4)/ ® I + (1 - C)|$probe)($probe|, (32) 

to see how necessary entanglement is for the QI—see 
Ref. [22] for a similar study in QPT and Ref. [30] for ex¬ 
perimental realization of Werner states. The parameter £ 
here adjusts the strength of quantum entanglement between 
paths siii and S 2 I 2 (see Fig. 3); W ext becomes entangled for 
0 ^ £ < 2/3 and separable for £ > 2/3 [31]. 

Running the quantum circuit (Fig. 2) by using the initial 
state (32) yields the following probabilities: 

Ps 1 /s 2 = (1/2) [l T (1 — £)T cos 7], (33) 

which (as expected) for the particular cases of £ = 0 and 
£ = 1 give (1/2) (1 =p T cos 7) and 1/2, respectively. This 
relation indicates that, in principle, any probe state with £ / 1 
can give an image. The quality of the image is quantified by 
its visibility, defined as (P max - Pmin)/(Pma X + P m in). For 
P si/S2 the visibility is equal to (1 — £)T, which increases 
by T and decreases by £. Thus one can obtain image even 
with non-entangled initial states. We remark that “quantum- 
mimetic” classical imaging has already been experimentally 
demonstrated; see, e.g., Ref. [32], where an image (of course 
with relatively smaller visibility) has been obtained by using 
classical light in an imaging scheme. 


V. SUMMARY 


We have demonstrated a detailed analogy between the 
quantum imaging method of Ref. [12] and ancilla-assisted 
process tomography. We have suggested a quantum cir¬ 
cuit corresponding to the experimental setup, which fully 
and identically reproduces the results of the imaging sce¬ 
nario and underlies how quantum imaging can be recast in 
quantum computation/information language. Through the 
tomography-imaging analogy we have analyzed the utility of 
non-entangled fields in imaging. In particular, we have ar¬ 
gued that it is possible to create an image if the maximally- 
entangled initial states of the AAPT version of the imaging 
setup is replaced by a separable Werner state. This implies 
that the quantum imaging scheme is not totally a result of 
quantum correlations in the form of entanglement. 

Our work can allow better understand behaviors of different 
identification or sensing schemes in the form a unified picture. 
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Appendix A: State preparation 

We start our circuit using four qubits in separated |0) modes as the input 

|*o) = |0) a i|0) < i|0) <2 |0) a2 , (Al) 

where “1” and “2” imply the source of photons (NL1 or NL2). A Hadamard gate followed by a CNOT play the role of the first 
beam splitter (BS1), which is responsible for the uncertainty about the path of the photon (Fig. 2). By applying on |<f>o), the 
Hadamard gate generates the following superposition: 

|$i> = (|0)si|0)ji|0)j 2 |0)s 2 + |0) s i|0)i 1 |l)j 2 |0) s2 )/v / 2 1 , 

which after the CNOT gate becomes 

l$2> = (|0) i i|l)ii|0) i2 |0) s2 + |0) sl |0)j 1 |l)i 2 |0) s2 )/v / 2 1 . (A2) 

At the next step, we apply a CNOT to simulate the first nonlinear crystal (NL1), 

1*3) = (|l)»i|l)ii|0) i2 |0) s2 + |0) sl |0) il |l) i2 |0) s2 )/V^, (A3) 

after which another CNOT is applied to complete the photon generation process, yielding Eq. (9). 


Appendix B: Object as a quantum channel 

In Ref. [12], a semi-transparent object is studied experimentally. We can attribute the following quantum channel (or opera¬ 
tion) to the object: 


£0: 


|0)d |0)d 

|l)d -► Te~* |l) d + VI^T^|1) W 


(Bl) 


Although this shows the correct action of the object, it does not represent an accurate mathematical description. The question 
is: what is a mathematically meaningful description of an object which transmits a photon with a factor ofTe 11 , reflects it with 
y/1 — T 2 , and does not give anything when there is no incident photon ? 

Here we argue that the object—as a scatterer of a field—can in fact be described as an amplitude-damping process , which can 
be modeled by a BS in the path of a photon [14, 33]. A BS couples the incident field mode d (represented by the associated 
annihilation operator ad) to another mode w (represented by the annihilation operator b w ), through the interaction H$b = 
z#(adb^ — a^b w ), with the corresponding unitary evolution Ubs = e lHsB [14, 34] (Fig. 6). 

We have Ubs &d U^ s = cos 0 e _2u;d ad + sin# e _2u;w b w , from whence 

f^Bs(|0)d|0)w) =|0)d|0)w, 
f^Bs(|l)d|0) w ) =^BS ( a d|0)d|0) w ) 

= cos 6 e^ d |l)d|0)w + sin# e* Ww |0) d |l)w- (B2) 


By tracing out over mode w we obtain the following Kraus operators for the operation of the object: 


K 0 = 


1 0 
0 T 


and K\ 


0 

0 0 ) ’ 


(B3) 


o 

(T, 7) 



FIG. 6. (color online). Object—simulated as a beam splitter—transmits the incident beam with probability cos 2 # and reflects it otherwise. 
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where \rj\ := cos 2 0 , and we have the trace-preserving condition K^Kq + K.\K\ = / [14, 33]. The object may also change the 
phase of the incident field; thus to be more general we assume y/rf = Te ~ 11 . 

The action of the Kraus representation £q [q\ = Yip KpQK^ on the input state 


f Qoo Qoi \ 
\ Qio Qn ) 


(B4) 


gives the final state as 


‘Eo[q] = 


Qoo 4- (1 — T 2 )qh Te n Qoi A 
Te ^ Ql o T 2 g n ) ' 


(B5) 


One can also find a more explicit representation for this process by expanding Kq and K\ in terms of the orthonormal basis 
operators {cro = /, o\ = cr x ,cr 2 = cr y ,cr 3 = cr z }/\/2, where cr a s are the Pauli matrices [35]. Specifically, defining the 

coefficients a 0a through Kp = E«=o a p<x°al'/%' yields < E 0 [q] = (1/2) T, a p(xo)a/3 VaQVp, where {xo) a p = E?=o a ia a ia 
can be considered a matrix representation of Eq in the orthonormal operator basis (l/\^){cr a }^ =0 . Thus we obtain 


|i + v^l 2 0 0 (i + VW-v^*) 

I V - v } 2 i Wi -1? 1 2 o 

0 i|\/T - v'l 2 IVI - V I 2 0 

V (i + VT*)(i - VT) o o Ii-VTI 2 


(B6) 


Note that although xo is a 4 x 4 positive-definite matrix with (1/2) ^ a p(xo)a(3&p&a = I (the trace-preserving condition), it 
only has 2 independent real parameters T and 7. 


Appendix C: Applying the object and MM operations 

Applying E^ [Eq. (10)] on the following state: 

Q2 =|^probe) (^probe I 

= |(| 1 >( 1 | ® | 1 )( 1 | ® | 0 )( 0 | 0 |0)(0| + 1 1)(o| 0 1 1)(0| ® |0)(1| ® |o)(1 1 

+ |0)(1| ® |0)(1| ® |1)(0| ® |1)(0| + |0)(0| ® |0)(0| ® |1)(1| ® |1)(1|) (Cl) 

yields 

g 3 =(X 0 £0 0 X® T)[q 2 ] 

= |(|1)<1| ® £b[|1)(1|] ® |0)(0| 0 |0)(0| +11)(0| 0 <Eo[ |1)(0|] 0 |0)(11 0 |0)(1| 

+ |0)(1| 0 £o[|0)(l|] 0 |1)(0| 0 |1)(0| + |0)(0| 0 £o(|0)(0|) 0 |1)(1| 0 |1)(1|) 

=l(r 2 |i)(i| 011)(1| 0 |o)(o| 0 |o)(o| + (1 - T 2 )|1)(1| 0 |o)(o| 0 |o)(o| 0 |o)(o| 

+ Te i7 |l)(0| 0 |1)(0| 0 |0)(1| 0 |0)(1| + Te _i7 |0)(l| 0 |0)(1| 0 |1)(0| 0 |1)(0| 

+ |0)(0| 0 |0)(0| 0 11)(11 0 |1)(1|). (C2) 

After applying MM [Eq. (1 1)], we obtain 

s =\ [r 2 |l)(l| 0 |S)(S| 0 1 0)(0| + (1 - T 2 )|l)(l| 0 |00)(00| 0 |0)(0| + Te i7 |l)(0| 0 |H)(H| 0 |0)(1| 

+ Te _ * 7 |0)(l| 0 |3)(3| 0 |1)(0| + 10)(0| 0 |3)(3| 0 |1)(1|], (C3) 

Now if we discard the idler photons, which means tracing out over i\i 2 , the state of the signal photons (sis 2 ) reduces to 

T = 1[|1)(1| 0 |0)(0| + Te* 7 |l)(0| 0 |0)(1| +Te- i7 |0)(l| 0 |1)(0| + |0)(0| 0 |1)(1|], (C4) 


Note that this is exactly the state qq [Eq. (7)] of the original experimental setup. 
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Appendix D: Bell-state measurements 

We recall the definition of the Bell states for the case of two qubits [14], 

l^ ±) = ^ (|00) ± |11 ) ) , (Dl) 

I**) =-L(|01> ± |10». (D2) 

It is straightforward to verify the following Schmidt decompositions for the Bell-state measurements: 

\$+)(<S>+\ = (II + XX-YY + ZZ)/4, (D3) 

\$-)($-\ = (II - XX+ YY + ZZ)/4, (D4) 

|^r+)($+| = (II + XX+ YY - ZZ)/A, (D5) 

|^r-)(^-| = (II - XX-YY - ZZ)/4. (D6) 


Appendix E: Measurements for the complete object tomography 

A phase shifter e %) has the following action on the X and Y Pauli matrices: 

X Z ^ = cos X + sin <p Y. 

Z^Y Z^ = — sin X + cos Y. 

As a result, one can see that adding Z^ before the BS and detectors in part B of Fig. 4 changes fWh/g 1° 
^h/g = (A ® I CNOT H 0 I CNCl) M h/g (cNOT H 0 I CNOT Zj, 0 /j 
( = \(Zt 0 I) [(II - ZZ) T (XX + YY)\ (Z\ 0 I) 

XY-YX\ 1 

Hence we see that this phase shifter allows the missing measurement XY — YX appear. 


(E1),(E2) \[2 \II - ZZ 


Vs 


COS ( 


XX+ YY 


— smc 


(El) 

(E2) 


(E3) 
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